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We report measurements of the charge-separated — > e +( -^ + v e (i> e ) and Z/-y* — ¥ e + e~ 
production cross sections at mid-rapidity in proton-proton collisions at ^/s = 500 GeV. These results 
are based on 13.2 pb -1 of data recorded in 2009 by the STAR detector at RHIC. Production cross 
sections for W bosons that decay via the ev channel were measured to be a(pp — > W + X) •BR(H /+ — > 
e+ v e ) = 117.3 ± 5.9(stat) ± 6.2(syst) ± 15.2(lumi) pb, and a(pp -> W~ X)-BR(W~ -> e" 9 e ) = 43.3 
± 4.6(stat) ± 3.4(syst) ± 5.6(lumi) pb. For Z/-y* production, a(pp Z/-/* X) ■ BR(Z/-y* -> e+e") 
= 7.7 ± 2.1(stat) to.g(syst) ± l.O(lumi) pb for di-lepton invariant masses m e + e - between 70 and 
110 GeV/c 2 . First measurements of the W cross section ratio, a(pp — > W + X)/a(pp — > W~ X), at ^/s 
= 500 GeV are also reported. Theoretical predictions, calculated using recent parton distribution 
functions, are found to agree with the measured cross sections. 
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I. INTRODUCTION 

Studies of inclusive W and Z/7* boson production 
in proton-proton collisions provide valuable information, 
both to test the Standard Model of particle physics and 
to advance our understanding of the proton's substruc- 
ture. Measurements of the production cross sections 
a(pp -> W+^X) ■ BR(VK+(- ) + v e (v e )) and 

a(pp -> Z/~f* X) ■ BR(Z/7* — > e+e~) can be compared 
to theoretical calculations that involve the weak cou- 
plings between intermediate vector bosons and quarks, 
and which must account for higher-order terms in per- 
turbative QCD. Such calculations also rely on models of 
the parton distribution functions (PDFs) for the quarks 
and, in pp collisions, for the antiquark 'sea.' 

Until recently, most measurements of W and Z/j* pro- 
duction in hadronic interactions have been confined to 
experiments using proton-antiproton collisions. First re- 
sults were obtained by the UAl [U \2\ and UA2 [3] g] 
collaborations at y/s = 630 GeV at the CERN SppS 
facility, followed by the CDF [B] and DO H] pp mea- 
surements at the Fermilab Tevatron, at y/s — 1.8 and 
1.96 TeV. It is only in the last few years that pp colliders 
have reached sufficient center of mass energies for com- 
parable studies, at y/s = 500 GeV by the STAR [5] and 
PHENIX [10] collaborations at the Relativistic Heavy Ion 
Collider (RHIC), and most recently by the LHC experi- 
ments ATLAS [TTJ and CMS QH [□] at V? = 7 TeV. 

RHIC is unique in its capability to collide high en- 
ergy polarized proton beams, and the observation of 
W production in these polarized proton collisions pro- 
vides a new means to explore the spin-flavor structure 
of proton sea quark distributions. First measurements 
of the parity-violating longitudinal single-spin asymme- 
try for W ± decay leptons have also been reported by 
the STAR [9] and PHENIX fTD] collaborations and are 
in good agreement with predictions from NLO and re- 
summed calculations [TH [15] . 

At hadron colliders, the leading process in pro- 
duction is u + d (d + u) fusion. This suggests that while 
the W + and W~ production cross sections should be 
close to equal in pp collisions, they can be expected to 
differ in pp measurements due to differences in the u and 
d quark and antiquark distributions within the proton. 
The PDFs that characterize the valence u and d quarks 
of the proton (or u and d in the antiproton) are well deter- 
mined from decades of high precision, deep-inelastic lep- 
ton scattering experiments (see, for example, Ref. |16j). 
Comparable distributions for the antiquarks within the 
proton sea, however, are much more weakly constrained. 
Interest in these poorly-known antiquark PDF's has also 
increased over the last few years, due to results from 
Drell-Yan experiments |17| 118] which find evidence for a 
much larger d/u flavor asymmetry in the nucleon than 
had been anticipated, especially at momentum fractions 
near and above x ~ 0.2. Detailed measurements of W 
and Z/7* production in proton-proton collisions will pro- 
vide new and complementary information about this fla- 



vor asymmetry in the sea, from different reactions and 
at very different momentum scales. 

This paper describes the first measurement of the W + , 
W~ , and Z l^i* boson production cross sections in proton- 
proton collisions at \fs = 500 GeV by the STAR collab- 
oration at RHIC. The cross sections are derived from 
studies of the charge-separated — > + v e {p^) 

and Z/7* — » e + e~ decay channels for outgoing leptons 
near mid-rapidity (\r] e \ < 1), and are based on 13.2 pb _1 
of data recorded during the 2009 run. In addition to 
the individual cross sections, a first measurement of the 
W + /W~ cross section ratio at y/s — 500 GeV is also 
presented. 

The paper is organized as follows. Section |H] provides 
a brief overview of the STAR detector, focusing on the 
subsystems used in this analysis. Section III describes 



the data and simulation samples analyzed, Sec. |IV| de- 
tails the extraction of the W and Z/-f* signal spectra, 
and Sec. [V] explains the estimation and subtraction of 
the background from the signal spectra. Finally, we dis- 
cuss the calcula tion of the W and Z/j* production cross 
sections in Sec. |VI| and the W + /W~ cross section ratio 
in Sec. VII and compare these results to several theo- 
retical calculations. Some of the data analysis methods 
employed here have been described briefly in Ref. [3] , and 
are discussed in more detail in this paper which incorpo- 
rates a slightly larger data sample as well as improved 
detector calibrations with respect to the previous publi- 
cation. 



II. THE STAR DETECTOR 

The STAR detector (Solenoidal Tracker at RHIC) [15] , 
shown schematically in Fig.[l] is a large acceptance, mul- 
tipurpose detector designed primarily for measurements 
of hadronic and electromagnetic particle production in 
high-energy heavy ion and polarized proton-proton col- 
lisions. STAR is comprised of many separate subsys- 
tems, each with specific capabilities; only those subsys- 
tems most relevant for the present analysis will be men- 
tioned below. 

The heart of STAR is a large Time Projection Cham- 
ber (TPC) [3U] which is situated within a highly uni- 
form, 0.5 T solenoidal magnetic field. The TPC provides 
charged particle tracking, particle identification (via ion- 
ization energy loss, dE/dx), and precision momentum 
measurements over the range |ry| < 1.3 and with full 27r 
azimuthal coverage. Although the pr resolution of the 
TPC deteriorates with increasing px, the spacial accu- 
racy of tracks reconstructed between the inner and outer 
radius of the TPC, located at 50 and 200 cm respectively, 
remains accurate up to ~l-2 mm in Cartesian space. In 
this analysis, TPC tracks were used in identifying the 
high-pT decay lepton (e ± ) candidates, determining can- 
didate charge signs, reducing contamination from the sig- 
nificant QCD background (see Sec. IV I, and reconstruct- 
ing the interaction vertex for the events of interest. 
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Surrounding the TPC radially is the Barrel Electro- 
magnetic Calorimeter (BEMC) [21 , a high granularity 
lead/scintillator-based sampling calorimeter. This detec- 
tor is used to measure the energy deposited by energetic 
photons and electrons with pseudorapidities \rj\ < 1.0 
over the full azimuth. The BEMC is segmented into 
4800 optically isolated projective towers, each of which 
subtends 0.05 rad in azimuth (</>) and 0.05 units in 77, 
and is roughly 20 radiation lengths deep. Based on 
cosmic ray and test beam data, the nominal energy 
resolution of the barrel calorimeter is calculated to be 
SE/E = 14%/ v /£(GeV) © 1.5% [21]. The BEMC was 
used to measure the e candidate energy, and to aid in 
background reduction. By identifying events with large, 
highly localized, electromagnetic energy deposition, the 
BEMC also provided our first-level trigger signal for lep- 
tonic W and Z decays. 

Located at one end of the STAR TPC, directly in 
front of the magnetic field return poletip, is the Endcap 
Electromagnetic Calorimeter (EEMC) [22], which pro- 
vides electromagnetic energy measurement over the range 
1.09 < 77 < 2 and 2ir in azimuth. The EEMC is simi- 
lar in design to the BEMC: a lead/scintillator sampling 
calorimeter, finely segmented in 77 and </> into 720 towers 
with projective geometries, though it is approximately 
3-4 radiation lengths thicker than the BEMC due to its 
more forward position. In the work presented here, the 
EEMC was used only as part of the background reduc- 
tion via isolation and missing energy conditions discussed 
in Sec. |lVl 

III. DATA AND SIMULATION SAMPLES 

Candidate events were selected online using a two-level 
trigger requirement in the BEMC. The hardware level-0 
trigger accepted events containing a tower with a trans- 
verse energy, Et, greater than 7.3 GeV. A dedicated 
software trigger algorithm then selected events by con- 
structing 2x2 clusters of towers, and requiring that at 
least one cluster consist of a seed tower with Et > 5 GeV 
and a cluster sum Et > 13 GeV. During the 2009 run 
1.2 x 10 6 events were recorded satisfying these trigger 
conditions. 

The integrated luminosity of the data sample was de- 
termined using the Vernier Scan technique [23] . The 
transverse widths (a x and a y ) of the beam overlap re- 
gion are determined by measuring the trigger rate as the 
beams are swept through each other in the transverse 
plane. The intensity of each beam is determined during 
a scan by the Wall Current Monitors (WCM) [23] . With 
the assumption of Gaussian beams, the instantaneous lu- 
minosity can be written as 

where f rev is the revolution frequency and K = N°"N^ 
is the product of the bunch intensities (iVj) of the two 




FIG. 1. (Color online) (a) W candidate event display embed- 
ded in a schematic of the STAR detector. Descriptions of the 
subsystems relevant for this analysis are given in Sec. [IT] (b) 
TPC p T and (c) BEMC and EEMC E T distributions in 77 and 
<j> for the same W candidate event as shown in (a). 



beams (a, b) summed over all bunches. The dedicated 
trigger used in the Vernier Scan, and also to monitor 
the luminosity in this analysis, is the level-0 hardware 
trigger, described above, with a coincidence away-side 
Et requirement imposed offline to reduce non-collision 
background. The cross section for this trigger can be 
written as a vcr — R™ x /£, where R™^ x is the maximum 
trigger rate while the beams are fully overlapping. The 
value measured for this work was <7 ver = 434 ± 8(stat) 
± 56(syst) nb. Figure [2] shows an example of the trigger 
rate as a function of the x and y beam displacements dur- 
ing one of the vernier scans, which was fit to extract the 
transverse beam widths and maximum trigger rate. The 
fit function used was a Gaussian in x and y combined with 
a constant term to account for remaining non-collision 
background. The largest contribution to the er vcr system- 
atic uncertainty was attributed to possible non-Gaussian 
components of the beam profile (10%), with smaller con- 
tributions coming from possible BEMC gain drift (5%), 
and uncertainties in the bunch intensity measurements 
(4%). This value for cr vor was used to normalize the to- 
tal number of events which satisfy this trigger condition, 
resulting in an integrated luminosity for the data sample 
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of L = J Cdt = 13.2 ± 0.2(stat) ± 1.7(syst) pb _1 . 




Beam displacement in y (mm) 

FIG. 2. (Color online) Trigger rate as a function of vernier 
scan beam displacement in the x and y directions. The trans- 
verse beam widths (a x and a y ) and maximum trigger rate 
(R™r X ) were extracted from the fit, which is superimposed. 



Monte-Carlo (MC) simulation samples were generated 
in order to determine detector efficiencies, estimate back- 
ground contributions from electroweak processes, and 
compare various predicted observables to data. Signal 
samples for both the W — > ev and Z/7* — > e + er chan- 
nels were generated, along with a,W^Tv sample which 
is an expected background in the W analysis due to the 
r's leptonic decay. All the samples were produced us- 
ing the PYTHIA 6.422 event generator and a GEANT 
[26] model of the STAR detector response. The same 
reconstruction and analysis algorithm was used for both 
the data and MC samples, and each MC sample was nor- 
malized to the integrated luminosity of the data unless 
otherwise stated. 

Due to the high luminosity of the pp collision environ- 
ment at y/s = 500 GeV at STAR, a significant number 
of pile-up tracks are present in the TPC at any given 
time. The pile-up tracks are the result of either another 
collision from the same bunch crossing as the triggered 
event, or a collision that occurred in an earlier or later 
bunch crossing. Note that the bunch crossing period at 
RHIC is about 107 ns, while it can take up to ~38 /is for 
track ionization to drift through the TPC. In the simu- 
lation, these pile-up tracks are accounted for by embed- 
ding the full GEANT detector response of the simulated 
event into a zero-bias triggered event before reconstruc- 
tion. The zero-bias events are selected randomly during 
nominal beam crossings at a rate of J$ 1Hz with no de- 
tector requirements, resulting in a good representation of 
the pile-up contained in the TPC for BEMC triggered 
collision events. 



IV. W AND Z/7* SIGNAL RECONSTRUCTION 

This section details the identification and reconstruc- 
tion of W and Z/7* candidate events, as well as the re- 
duction of the large QCD background. This reduction 
is achieved through a number of cuts designed to take 
advantage of the kinematical and topological differences 
between electroweak and QCD processes. "Z/7*" will be 
used interchangeably with "Z" for the remainder of this 
paper. 

Candidate events were selected from the sample of 
BEMC triggered events described in Sec. |III| by requir- 
ing a reconstructed primary vertex. A primary vertex is 
one reconstructed from either a single TPC track with 
Pt > 10 GeV/c or multiple tracks originating from the 
same location along the beamline. Each track considered 
in vertex reconstruction is assigned an increased weight 
if it either points to a region of energy deposition in the 
calorimeters, or if it uses hit points from both sides of 
the TPC central membrane. Tracks satisfying either of 
these two conditions are likely to be from the triggered 
collision, therefore weighting these tracks more heavily in 
vertex reconstruction strongly reduces the contamination 
from pile-up tracks. The distribution of primary vertices 
along the beam direction is approximately Gaussian with 
an RMS width of 52 cm. Events of interest were required 
to have a | Vertex I < 100 cm, where z ver t GX is the distance 
along the beam direction of the primary vertex from the 
nominal collision point at the center of the STAR inter- 
action region. 

A. Identification of High-iSr Isolated Electrons and 
Positrons 

A candidate electron or positron track is defined to be 
a TPC track with pt > 10 GeV/c that is associated with 
a primary vertex satisfying the criteria described above. 
Candidate tracks were also required to have: 

• a minimum of 15 TPC points, 

• more than 51% of the maximum number of TPC 
points allowed, 

• a first TPC point with radius less than 90 cm, 

• a last TPC point with radius greater than 160 cm. 

These requirements help to ensure that the track and its 
charge sign are well reconstructed, as well as reject pile- 
up tracks which may be mistakenly associated with a 
primary vertex. Candidate TPC tracks are extrapolated 
to the BEMC to determine to which tower the track 
points, then the four possible 2x2 BEMC tower clusters 
containing the tower pointed to by the track are formed. 
The 2x2 cluster with the largest summed transverse en- 
ergy, Ej*, is assigned to the candidate. The candidate 
Efp is required to be greater than 15 GeV to be safely 
above the trigger turn-on region. Also, the two dimen- 
sional distance between the energy log-weighted centroid 



of the tower cluster position and the extrapolated TPC 
track position, |Ar|, is required to be less than 7 cm, to 
reject candidates where the BEMC cluster may not have 
originated from the particle which produced the high px 
TPC track. 

Electrons and positrons from W and Z decays should 
be well isolated from other particles in rj — <fi space; thus, 
in the next stage of the candidate selection process two 
isolation criteria are applied. The first isolation cut was 
made by summing the Et in the 4x4 BEMC tower clus- 
ter which surrounds the e ± candidate cluster, E T x4 , and 
requiring E T /E T xi > 0.95. The other isolation require- 
ment is imposed to reduce jet-like events. The quantity 
e£ r<0 - 7 is defined as the sum of all BEMC and EEMC 
tower Et and TPC track pt within a cone radius of 
AR = yj At] 2 + Acj) 2 < 0.7 around the candidate track, 
and the ratio E T /E^ K<0 ' 7 is required to be greater than 
0.88. The e~^ candidate track is excluded from the sum 
of TPC track pt to avoid double-counting the candidate 
energy in the E^ K<0 ' 7 sum. Figure shows the isola- 
tion ratios described above for both data and W — » ev 
MC. The placements of the cuts, shown by the dashed 
lines, were chosen to retain a large fraction of the signal, 
while significantly reducing the QCD background. Note 
that differences between the isolation ratios in Fig. [3] of 
this paper and Fig. 1 of Ref. [S] are expected due to 
differences in the data samples used and improved cali- 
brations. Also, the order of the E T /E T xA and candidate 
track-cluster matching |Ar| cuts were inverted in Rcf. [9 
with respect to the ordering described in this section. 
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FIG. 3. (Color online) Distributions of the isolation ratios 
E^/E^ xi (left) and E^/E^ R<0 - 7 (right) used in e ± candi- 
date selection. W — > ev MC shape distributions (arbitrary 
normalization) are shown as filled histograms for comparison 
with the data distributions. The vertical dashed lines indicate 
the placements of the cuts on these isolation ratios. 



B. W Candidate Event Selection 

The selection of W — > ev candidate events is based 
on differences in the event topology between leptonic W 
decays and the QCD background or Z events. W — > ev 
events contain a nearly isolated e ± with a neutrino close 
to opposite in azimuth. Electrons and positrons emitted 
near mid-rapidity from W decay are characterized by a 
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large E T that is peaked near half the W mass (~40 GeV) 
with a distribution referred to as a Jacobian peak. There 
is also a large missing transverse energy in W — » ev 
events opposite, in azimuth, to the e due to the unde- 
tected neutrino. As a result, there is a large imbalance in 
the vector pt sum of all reconstructed final state objects 
for W events. In contrast, Z — > e + e~ events and QCD 
hard-scattering events, such as di-jets, are characterized 
by a small magnitude of this vector pt sum imbalance. 

In order to enforce this missing energy requirement, we 
define the pt balance vector: 

P* al =P T e + E PT"' 8 (2) 
AR>0.7 

where p^ is the e ± candidate pt vector, which is com- 
posed of a momentum direction and a magnitude deter- 
mined by the candidate TPC track and BEMC cluster, 
respectively. The second term on the right of Eq.[2]is the 
sum of the pt vectors for all reconstructed jets whose 
thrust axes are outside the cone radius of AR = 0.7 
around the candidate. Jets are reconstructed using a 
standard mid-point cone algorithm used in STAR jet 
measurements j27j based on the tracks from the TPC 
and tower energies in the BEMC and EEMC. A scalar 
signed Pp-balance variable is then formed, defined as 

signed P T -balance = sign (p T e • p^ al ) \p^ al \ . (3) 

This quantity is required to be larger than 15 GeV/c as 
indicated by the dashed line in Fig. |4j Also in Fig. [4] 
one can see that in the W — > ev MC sample, the signed 
Pr-balance variable and E T are very well correlated, as 
contributions to the Pt"- 1 vector from reconstructed jets 
outside the cone of AR = 0.7 are generally small. The 
data show a similar correlation at high E T , where the 
distribution is dominated by W — > ev events. At low 
E T where contributions from QCD background events 
are larger, more events have a small value for the signed 
Py-balance variable, as expected. 




20 40 60 20 40 60 



(GeV) E?(GeV) 

FIG. 4. (Color online) Correlation of the signed Pr-balance 
variable and for data (left) and W — > ev MC (right). 

Background events from Z — > e + e~ decays are further 
suppressed by rejecting events with an additional e-like 



7 



2x2 cluster in a reconstructed jet where E^ x2 > /2 
and the invariant mass of the two e ± -like clusters is 
within the range of 70 to 140 GeV/c 2 . This reduces 
Z — > e + e~ contamination in both the W signal spectra 
and in the spectra that will be used for the data-driven 
QCD background, described in Sec. |V A| 

The reduction in the W candidate yield after each of 
the selection criteria is shown in Fig. [5] Initially, when 
only a candidate TPC track and BEMC cluster have 
been reconstructed, the distribution (solid line) is domi- 
nated by QCD background, which is exponentially falling 
with Ej,, and there is no evidence of the Jacobian peak. 
However, once the selection, isolation and signed Pt~ 
balance cuts are applied, a W signal can be seen above 
the background at E^ ~ Mw/2- 



Candidate track and BEMC 

cluster reconstructed 
E^> 15GeVand |Ar| < 7 cm 

E^/E 4 .* 4 > 0.95 

E?/E* R<0 ' 7 > 0.88 

Signed P T - balance > 15 GeV/c 



Z Candidate Event Selection 




40 

E*(GeV) 

FIG. 5. (Color online) Distributions of E^ for W candidate 
events after sequentially applying the selection criteria de- 
scribed in Sees. HVAl and HTBl 

The charge sign of the candidate is determined by 
the direction of curvature of the TPC track in the STAR 
magnetic field, while the magnitude of the track curva- 
ture provides a measure of l/pr- Figure [6] shows the 
product of the reconstructed charge sign and l/pr for 
the lepton candidates that satisfy all the cuts described 
above with Ej, > 25 GeV. Two well-separated regions 
are seen for the positive and negative charges, cleanly 
distinguishing between the e + and e~ candidates. 




0.02 0.04 0.06 
Lepton charge sign x 1/p (c/GeV) 



FIG. 6. (Color online) Distribution of the product of the TPC 
reconstructed charge sign and 1/pt for candidates satisfying 
all the W signal selection criteria and E^ > 25 GeV. 



Using the isolated e ± sample found in Sec. |IV A 
Z — > e + e~ events were selected by requiring a pair of 
isolated e ± candidates with opposite charge signs. The 
invariant mass of each e + e~ pair was reconstructed, and 
the resulting mass distributions are shown in F ig. [7] a fter 
each of the selection criteria described in Sec. IIVAI has 
been satisfied for both the e + and e~ candidates. After 



Candidate track and BEMC 

cluster reconstructed 
E!f > 15 GeV and |Ar| < 7 cm 
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E^/E 4 .* 4 > 0.95 




ny e - (GeV/c 2 ) 



FIG. 7. (Color online) Distributions of the invariant mass of 
pairs of oppositely charged e candidates after sequentially 
applying the selection criteria described in Sec. |IV A] to both 

e 



all selection cuts are applied, there is a signal near the 
invariant mass of the Z and a small signal at lower in- 
variant mass. This is consistent with the expectations 
from the Z/7* e+e~ MC, as shown in Fig. [Sj 




m e , e (GeV/c 2 ) 



FIG 

Z/j* -> e+e 
described in Sec 



(Color online) Distributions of the invariant mass of 
ca ndidat e events satisfying all selection criteria 
IV C| The Z/7* e + e~ MC distribution 
(dashed line) is shown for comparison. Note the larger bin 
widths relative to Fig. [7] 



V. BACKGROUND ESTIMATION 
A. W Background Estimation 

There are a number of background processes that can 
contribute to the W — > ev candidate yield. Other elec- 
troweak processes also yield isolated electrons that can 
be misidentified as W — > ev events, and QCD jets can 
fragment in such a way that they satisfy all the W signal 
requirements. This section describes how the contribu- 
tions of these background processes to the W candidate 
yield are estimated. 

The electroweak background processes considered in 
this analysis are W — > tv and Z — > e + e~ . Their contri- 
butions to the W — > ev signal yield were estimated using 
the MC samples described in Sec. |III| W — >■ tv events, 
where the r decays leptonically {i.e. r — > evv), contain 
an isolated e with a large missing energy opposite in 
azimuth, similar to the W — > ev signal. However, the 
e which comes from the r decay must share the energy 
of the r with the two secondary neutrinos, and thus it 
has a much lower Ej, on average than those e which 
come directly from a W decay. Therefore, the W — » tv 
background contributions are largest at low Ej,, as can 
be seen in Fig. [9J Z — > e + e~ events can contaminate 
the W signal when one of the decay leptons escapes de- 
tection, either from a detector inefficiency or by emission 
into an uninstrumented region of phase space. Unlike the 
other background sources described in this section, the 
Z — > e + e~ background yield is approximately constant 
in Ej,, resulting in a significant contribution to the to- 
tal background, even though the cross section is small 
compared to other processes. Table |T] lists each of the 
background processes and its estimated contribution to 
the W yield for candidates with Ej, > 25 GeV. The un- 
certainties for these electroweak background components 
are due to the statistical uncertainty of the MC calcula- 
tion and the uncertainty in the normalization of the MC 
samples to the integrated luminosity of the data. 

The STAR detector has only one EEMC, resulting 
in missing calorimetry acceptance for the pseudorapid- 
ity region —2 < r/ < —1.09 compared to the positive 
pseudorapidity portion of the detector. If the isolation 
cone of AR < 0.7 around an e ± candidate overlaps with 
this missing acceptance, or a jet opposite in azimuth 
of an e candidate falls within this acceptance, back- 
ground QCD events may satisfy all the W — > ev selec- 
tion requirements. This contamination of the W yield, 
referred to as the 'second EEMC background, was de- 
termined by repeating the W signal selection a second 
time, with the EEMC towers excluded from the isola- 
tion ratio, Ej,/e£ r<0 ' 7 , and the reconstruction of jets 
summed in the prf al vector. The events which satisfy the 
requirements of this second pass analysis (without the 
EEMC), but fail the nominal requirements described in 
Sees. IIV Al and IIVBI arc a direct measure of the back- 
ground rejected by the EEMC. Moreover, these events 
also estimate the amount of background that would have 



been rejected by a second EEMC. 

While this sample of second EEMC background is ex- 
pected to be predominantly the result of QCD processes, 
it does contain a small amount of Z — > e + e~ contami- 
nation as well. Because background from the Z — > e + e~ 
process was already taken into account separately, the 
Z — > e + e~ MC sample was used to remove any contam- 
ination from Z — > e + e~ processes in the second EEMC 
background distribution, to avoid double-counting. The 
uncertainty on the second EEMC background is the sta- 
tistical uncertainty of the events vetoed by the EEMC 
and the systematic uncertainty in the normalization of 
Z — > e + e~ contamination which was subtracted using 
the Z -> e+e- MC. 

The remaining contribution to the background is pre- 
dominantly from QCD 2 — > 2 processes in which one jet 
fragments such that it satisfies our e candidate require- 
ments, while all other jets escape detection outside the 
\rj\ < 2 acceptance. This component of the background 
was estimated using a data-driven QCD background dis- 
tribution as a function of Ej,, which is obtained by se- 
lecting events which satisfy all the isolated e candidate 
criteria, but have a signed Py-balance < 15 GeV/c. Sim- 
ilar to the way the second EEMC background was cor- 
rected, contributions to the data-driven background dis- 
tribution from the Z — > e + e~ process were removed using 
the Z — > e + e~ MC sample, to avoid double-counting the 
Z — > e + e~ background. 

The data-driven QCD background distribution was 
then normalized to the remaining W — > ev candidate 
signal distribution after the W — > tv, Z —> e + e~, 
and second EEMC background components had been 
removed. The normalization was determined over the 
range 15 < Ej. < 19 GeV, and accounts for the possi- 
bility of true W signal events in this region using the 
W — > ev MC. The systematic uncertainty of this data- 
driven QCD background contribution was estimated by 
varying the data-driven background distribution and the 
E^ region over which the distribution was normalized. 
Twenty different background distributions were obtained 
by varying the cut on the signed Py-balance variable 
from 5 to 25 GeV/c in steps of 1 GeV/c. The twenty 
background distributions were then fit to the signal, as 
described above, using three different normalization re- 
gions (15 < Ej, < 17, 19, and 21 GeV), resulting in sixty 
different normalized background distributions. The sys- 
tematic uncertainty in each Ej, bin was taken to be the 
largest deviation among these sixty distributions from 
the nominal value. 

The charge-separated Efp distributions of W ± — > e + 
v e candidates satisfying all the selection criteria described 
in Sees. |IV A| and |IV B| are shown in Fig. [9| Also shown 
here are the contributions from the different backgrounds 
discussed in this section and the W — > ev signal MC dis- 
tribution. A x 2 test of homogeneity comparing the data 
and the sum of background components and W — > ev 
signal MC (dashed line) Ej, spectra results in a \ 2 value 
of 9.5 and 6.9 for the W + and W~, respectively. For 12 
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Positron |n I < 1 

STAR 2009 Data 

— W-)ev MC 
Data-driven QCD 




Electron |r|J < 1 
Second EEMC 

■ W — > x v MC 

■ Z -> ee MC 



30 40 
E?(GeV) 

FIG. 9. (Color online) E% distribution of W + (top) and 
W~ (bottom) candidate events, background components, and 
W — > ev MC signal for comparison. Note the factor of two 
difference in the vertical scales. 





W+ -¥ e+ v a 


W — > e u e 


W — > TV 

Z -> e+e" 
Second EEMC 
Data-driven QCD 


13.4 ± 1.7 ± 3.2 
7.3 ± 0.4 ± 1.7 
9.1 ± 3.0 ± 0.5 
7.0 ± 0.6 tH 


3.3 ± 0.8 ± 0.8 
7.3 ± 0.4 ± 1.7 
9.2 ± 3.0 ± 0.4 
5.8 ± 0.5 t 2 { 6 2 


Total 


36.6 ± 3.5 ttt 


25.8 ± 3.2 ±f;f 



TABLE I. Summary of background event contributions to the 
W — > ev yield and their uncertainties for candidates with 
E e T > 25 GeV and \rj e \ < 1. 
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FIG. 10. (Color online) Lepton pseudorapidity distribution 
of W + (left) and W~ (right) candidate events, background 
components, and W — > ev MC signal for comparison. 



in Sec. Mil Within the defined mass window to be used 
for the cross section (70 < m e + e - < 110 GeV/c 2 ), the 
background contributions were determined to be 0.1 Iq'i 
events from W — > ev and negligible from the other Z 
decay channels. The W — > ev background uncertainty 
was estimated using the 68% C.L. interval of the unified 
statistical approach described in Ref. [55] , 

An accurate data-driven estimate of the QCD back- 
ground is difficult to obtain for the Z signal due to the 
limited statistics of the data set. One method for esti- 
mating the background is to determine the number of e ± 
pairs that satisfy all the Z — > e + er signal criteria other 
than the opposite charge-sign requirement. However, no 
same charge-sign pairs were observed in the data, there- 
fore, the QCD background was found to be consistent 
with zero. An upper bound on the QCD background 
systematic uncertainty was estimated to be 1.3 events 
using a 68% C.L. interval [IB]- 



VI. THE W AND Z CROSS SECTIONS 



degrees of freedom this results in a 66% and 86% proba- 
bility, respectively, to obtain a larger v 2 . This indicates 
a good agreement between data and MC and further val- 
idates the procedure used in the background estimation 
described in this section. The pseudorapidity distri- 
butions are shown in Fig. |10[ where the background con- 
tributions were found independently for each r/ e bin using 
the methods described above. Again, good agreement is 
found between the data and the sum of the W — > ev 
signal MC and background components. 



The W and Z production cross sections were measured 
from the sample of events which satisfy the fiducial and 
kinematic requirements of this analysis. As stated previ- 
ously, only e ± candidates at mid-rapidity (|rj e | < 1) were 
considered in this analysis. Candidates for the W anal- 
ysis must have Ej, > 25 GeV, and for the Z analysis we 
required that both e + and e~ have Ej, > 15 GeV and 
70 < m e + e - < 110 GeV/c 2 . The cross sections measured 
within these constraints are defined as the fiducial cross 
sections, and can be written as: 



B. Z Background Estimation 

The background for the Z — > e + e~ signal is expected 
to be very small due to the coincidence requirement of 
a pair of oppositely charged, high Et, isolated e + and 
e~. Background contributions from electroweak pro- 
cesses were estimated using the MC samples described 



fid 
T W(Z) 



BR(W(Z) -> ev{ee)) 



Ajobs 



where 



Ajbkgd 
ly W(Z) 



T Aot 
^ ' e W(Z) 



(4) 



^w(z) ^ s t ne number of observed W(Z) candidates 
within the defined kinematic acceptance, which sat- 



isfy all the selection criteria described in Sec. IV 
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N w ? z s is the total number of W(Z) background 
events within the defined kinematic acceptance de- 
scribed in Sec. |Vl 



W(Z_ 



is the total efficiency correction described in 



Sec. VIA below 



• and L is the integrated luminosity of the data set 
discussed in Sec. Mil 

To determine the total production cross sections, it is 
necessary to apply acceptance correction factors, Aw^z), 
to the fiducial cross sections defined above, to account 
for the fiducial and kinematic constraints imposed in the 
analysis. The total production cross sections are then 
defined via the relations 



(7^* • BR(IT -> 



fid 



BR(FT 



A 



w 



BR{Z 



BR(Z 



(5) 



(0) 



The determination of the acceptance corrections neces- 
sary to extract the total production cross sections is dis- 
cussed in Sec. I VI Cl 



A. The Efficiency Correction Factors 

The efficiency corrections were obtained using the 
W —r ev and Z — f e + e~ pythia MC samples described 
in Sec. |III| Only the subset of events from the MC sam- 
ples which satisfy the acceptance conditions for the fidu- 
cial cross sections were used in the efficiency calculations, 
as the acceptance correction is accounted for separately 
in the definition of the total cross section. 

The total efficiency can be factorized into four condi- 
tional efficiency terms, written as: 



trig 

— C W(Z) 



Ark 



algo 



-W(Z) C W{Z) C W(Z)- 



(7) 



The values for each of the terms in Eq. [7] are listed in 
Table [TT1 along with their uncertainties, for the W + , W~ , 
and Z signals. The remainder of this section describes 
how those values were obtained. 

The trigger efficiency, e'™ 9 , is the fraction of MC signal 
events which satisfy the online trigger condition defined 
in Sec. |III| This was determined by emulating the trigger 
condition used online in the MC. Due to the relatively 
wide z vor tox distribution of our data sample, some candi- 
dates may satisfy the \~q e \ < 1 kinematic condition at the 
MC generator level, but will fall outside the acceptance 
of the BEMC. This was observed in the W analysis as an 
-Ef.-dependent trigger efficiency due to the correlation of 
the Ej* and r\ e of the decay e ± . An independent trigger 
efficiency correction was therefore used in the computa- 
tion of the W ± cross sections. This effect also leads to 



a notably smaller average W~ trigger efficiency relative 
to W + , as the i] e distribution is expected to be peaked 
more strongly at zero for the W + candidates than W~ , 
which is consistent with Fig. [TU] To estimate the uncer- 
tainty on e*"s, the BEMC energy scale was varied by its 
uncertainty of ±3.6%. Because the offline kinematic re- 
quirement of Ej, > 25 GeV was significantly larger than 
the trigger threshold of 13 GeV, for this analysis we ob- 
served only small variations in the trigger efficiency due 
to the uncertainty of the BEMC energy calibration. 

The vertex efficiency, e vert , is defined as the frac- 
tion of events satisfying the trigger which contain a re- 
constructed primary vertex within the fiducial cut of 
l^vertex! < 100 cm, as described in Sec. |lVl The track- 



ing efficiencies for the W and Z decay e ± s are defined as 
follows. For W events with a reconstructed primary ver- 
tex, e^ fe is the efficiency for reconstructing a single TPC 



track which satisfies the track requirements in Sec. |IV A 

c trk 
Z 



however for Z — > e + e~ events the tracking efficiency, e\ 
is the efficiency for reconstructing two TPC tracks sat- 
isfying those conditions. In comparing the reconstructed 
TPC track 1/pr distributions between data and MC, a 
slightly worse resolution was seen in the data. This was 
accounted for by re-weighting the MC distributions to 
match the data. The uncertainty on the tracking effi- 
ciency was estimated from the error in this re-weighting 
resulting from the limited statistics of the data distribu- 
tion. 

Finally, the algorithm efficiency, e al9 ° ', is the fraction of 
events with one (two) reconstructed e^ candidate TPC 
tracks, which satisfy the remaining W (Z) selection cri- 
teria. As discussed in Sec. |IV[ these remaining selection 
criteria include reconstruction of BEMC clusters, match- 
ing extrapolated track and cluster positions, isolation re- 
quirements, and finally the signed P-p-balance and pair of 
opposite charge-sign candidate requirements for W and 
Z events, respectively. A weak Ej, dependence was ob- 
served in the algorithm efficiency for the W — > ev MC 
due mainly to the efficiency of the i?f./E^ R<0 ' 7 isolation 
cut being reduced at low Ej,. Thus, an i?|,-dependent 
algorithm efficiency correction was used in the computa- 
tion of the cross sections. The uncertainty on e al9 ° 
was determined by varying the BEMC scale uncertainty, 
as was done for the trigger efficiency. 





W + -> e + u e 


W — ¥ e v e 


Z -¥ e+e" 


triq 

trk 
alqo 


0.857 ± 0.007 
0.881 ± 0.005 
0.741 ± 0.030 
0.892 ± 0.024 


0.825 ± 0.007 
0.886 ± 0.006 
0.748 ± 0.031 
0.892 ± 0.024 


0.968 ± 0.006 
0.938 ± 0.006 
0.511 ± 0.032 
0.730 ± 0.024 


tot 


0.498 ± 0.026 


0.488 ± 0.026 


0.338 ± 0.024 



TABLE II. Summary of conditional efficiency correction fac- 
tors included in Eq. [7] The average values for the trigger and 
algorithm efficiencies for the W analysis are given here, how- 
ever an _E^-dependent correction was used for the measured 
cross section, as described in the text. 
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B. The Measured Fiducial Cross Sections 

The fiducial cross sections are calculated according to 
Eq . [4j and the measured values are summarized in Tables 
IIl|and [rv] for W ± and Z respectively. The dominant un- 
certainty for both the W + and W~ cross sections comes 
from the systematic uncertainty in the measured lumi- 
nosity of the data sample. The Z cross section measure- 
ment, however, is currently dominated by the statistical 
uncertainty. 





W+ -> 


e + u e 




W~ -> 


e v e 






value 


stat 


syst 


lumi 


value 


stat 


syst 


lumi 


N obs 


496 


22.3 






148 


12.2 






jybkgd 


36.6 


3.5 


+5.4 
-5.2 




25.8 


3.2 


+3.6 
-2.8 




£ tot 


0.498 


0.006 


0.025 




0.488 


0.007 


0.025 




L (pb" 1 ) 


13.2 


0.2 




1.7 


13.2 


0.2 




1.7 


G fli ( P b) 


70.0 


3.5 


3.5 


9.1 


19.2 


2.1 


1.1 


2.5 



TABLE III. Summary of input and measured values for the 
W — > ev fiducial cross sections, with their statistical, system- 
atic, and luminosity uncertainties. As noted in the text, an 
Z?y-dependent efficiency correction factor is used for the cross 
section measurement, and only the average value is shown 
here. 





Z -> e+e" 




value 


stat 


syst 


lumi 


N° ba 


13 


3.6 






jy-bkgd 


0.1 


0.1 


+ 1.3 
-0.0 




£ tot 


0.338 


0.012 


0.021 




L (pb" 1 ) 


13.2 


0.2 




1.7 


<j fid (pb) 


2.9 


0.8 


+0.2 
-0.3 


0.4 



TABLE IV. Summary of input and measured values for the 
Z — > e + e~ fiducial cross section, with their statistical, sys- 
tematic, and luminosity uncertainties. 



C. The Acceptance Correction Factors 



As stated previously, to determine the total cross sec- 



tions, acceptance correction factors A 



W(Z) 



must be used 



to account for the fiducial and kinematic acceptance re- 
quirements of the analysis, which are defined at the be- 



ginning of Sec. VI Aw(z) were calculated using the fewz 
program which provides cross section calculations 
for W and Z boson production up to NNLO in pQCD. 
Table M lists the values of the acceptance factors using 
the MSTW 2008 [3U] and CTEQ 6.6 [3T] parton distri- 
bution function sets. The nominal values for the accep- 
tance corrections, used in the total cross section measure- 
ments, were taken from the next-to-leading order (NLO) 
calculation using the MSTW08 PDF set. Theoretical un- 



certainties in the calculation of these factors arise from 
several sources, including differences between PDF sets, 
uncertainties within a PDF set, and uncertainties in the 
modeling of the production process. 





A w + 


A w - 


A z 


LO MSTW08 


0.591 


0.444 


0.377 


NLO MSTW08 


0.597 


0.444 


0.378 


NNLO MSTW08 


0.603 


0.435 


0.385 


NLO CTEQ 6.6 


0.592 


0.430 


0.370 



TABLE V. Summary of acceptance values calculated with 
the fewz program. The NLO MSTW08 values are used for 
the total cross section calculations in Sec. I VI Dl 

The uncertainty due to differences between PDF sets 
was taken to be the difference between the CTEQ 6.6 and 
MSTW08 acceptance values at NLO. Both groups pro- 
vide error eigenvector PDF sets which were used to esti- 
mate the acceptance uncertainty, at the 90% confidence 
level, within each set. The average of the CTEQ 6.6 and 
MSTW08 error eigenvector uncertainty was taken to be 
the uncertainty due to the PDF itself. Finally, the uncer- 
tainty in the modeling of the production process was esti- 
mated by comparing the acceptance values from calcula- 
tions with different orders of QCD corrections, using the 
MSTW08 PDF set. The maximum difference from the 
nominal value (NLO MSTW08) was taken as this final 
uncertainty contribution. Table [VI] summarizes the con- 
tributions to the uncertainties in the acceptance values. 
The individual contributions were added in quadrature to 
determine the total uncertainty for each acceptance fac- 
tor. The A w - uncertainties are significantly larger than 
those for A w + , driven primarily by the PDF-related er- 
rors. This is expected, due to the larger uncertainties in 
the u and d quark PDFs with respect to those of the d 
and u quarks. 





5A w+ (%) 


SA w -(%) 


SA Z {%) 


Difference between PDFs 


1.0 


3.2 


2.1 


MSTW08 NLO error PDFs 


0.9 


2.7 


1.2 


CTEQ 6.6 NLO error PDFs 


0.9 


4.5 


1.8 


Calculation Order 


1.0 


2.0 


1.9 


Total 


1.7 


5.2 


3.2 



TABLE VI. Summary of the relative uncertainties in the 
acceptance correction factors, A W ( Z ), as computed by the 
fewz program. 



D. The Measured Total Cross Sections 

The total cross sections are calculated according to 
Eqs. [5] and [6j by dividing the measured fiducial cross 
sections by the acceptance correction factors determined 
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in the previous section. The results for pp —> W to- 
tal production cross sections at y/s — 500 GeV are the 
following: 

■ BR(W+ -> e+ v e ) = 117.3 ± 5.9(stat) 
± 6.2(syst) ± 15.2(lumi) pb, 

a w L • BR(W~ -> e- P e ) = 43.3 ± 4.6(stat) 
± 3.4(syst) ± 5.6(lumi) pb. 

The result for the pp — > Z/j* total production cross sec- 
tion at yfs — 500 GeV in the invariant mass range of 
70 < m e+e - < 110 GeV/c 2 is 

a^f T ■ BR(Z/ 7 * e+e") = 7.7 ± 2.1(stat) 
toj(syst) ± l.O(lumi) pb. 

Figure [TT] shows the measured total cross sections, mul- 
tiplied by the respective branching ratios, in comparison 
with the theoretical predictions at NLO from the fewz 
program using the MSTW08 PDF set. Measurements 
from other experiments at the SppS, Tevatron, RHIC, 
and LHC are also shown as a function of y/s for compar- 
ison. 



the rhicbos calculation. Variations in the strong cou- 
pling constant, a s , from the associated error PDF sets 
were considered as well, but the uncertainties were found 
to be negligible compared to the uncertainties from the 
PDFs. The theoretical predictions agree well with the 
measured cross sections within the theoretical and exper- 
imental uncertainties. Interestingly, differences between 
the MSTW08 and CTEQ 6.6 PDF sets result in signifi- 
cant differences in the predicted cross sections at NLO. 





<+(pb) 


<C-(pb) 


4°' • (Pb) 


NLO MSTW08 
NNLO MSTW08 


132.4 ± 9.0 
136.7 ± 9.5 


45.7 ± 3.6 
48.1 ± 3.0 


10.8 ± 0.8 
11.2 ± 0.8 


NLO CTEQ 6.6 


121.8 ± 8.8 


41.1 ± 4.3 


9.8 ± 0.8 


Ressum. CTEQ 6.6 


121.1 


39.9 





TABLE VII. Summary of total cross section (times branching 
ratio) theoretical predictions at i/s — 500 GeV calculated 
with the fewz and RHICBOS programs. The Z/j* values are 
defined within the invariant mass range of 70 < m e + e - < 
110 GeV/c 2 . 
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FIG. 11. (Color online) Measurements of W and Z total 
cross sections times branching ratio versus center-of-mass en- 
ergy. For the W cross sections in pp collisions, the closed 
symbols represent W + and the open symbols represent W~ . 
The theory curves are from the fewz program at NLO using 
the MSTW08 PDF set. 

Theoretical predictions for the production cross sec- 
tions computed by the fewz [55] and fully r esum mcd 
RHICBOS [15] calculations are shown in Tabic VII The 
theoretical uncertainties were determined for the fewz 
predictions using the 90% confidence level error eigenvec- 
tor PDF sets; error eigenvector sets are not provided for 



VII. THE W CROSS SECTION RATIO 



The W cross section ratio is defined as 



Rw = 



fid 

a w+ 


N obs 


nj-bkgd 

iy w+ 


t w- 


fid 

°w- 


ATobs 

w- 


Ajbkgd 

w- 





(8) 



If the small contributions from strange quarks are ne- 
glected, this ratio should be equal to [3"2"] 



Rw 



u(xi)d(x 2 ) + d{xx)u(x2) 
u(xi)d(x 2 ) + d(xi)u(x 2 ) ' 



(9) 



Measurements of the cross section ratio should therefore 
be sensitive to the flavor asymmetry of the antiquark sea 
in the Bjorken-x range 0.1 < x 0.3 probed at RHIC. 
Drell-Yan experiments [TTl [18] have measured a large 
asymmetry in this x range, and precision measurements 
of Rw at RHIC can provide independent constraints on 
the flavor asymmetry which are free from the assump- 
tion of charge symmetry required in Drell-Yan. Mea- 
surements of the lepton charge asymmetry at the LHC 
[331 134j provide similar constraints on the quark and an- 
tiquark PDFs, though at significantly lower x due to the 
much higher energy of the collisions. 

The W cross section ratio was measured in two \rj e \ re- 
gions, as this coarsely constrains the x of the partons in- 
volved in the W production. In each \rj e \ bin, the fiducial 
cross sections were computed using the same procedures 
described in Sec. |VI[ where the background and efficien- 
cies were separately calculated for each charge and \rj e \ 
bin. The luminosity, and its sizable uncertainty, cancel in 
the cross section ratio, significantly reducing the system- 
atic uncertainty, with respect to the W + and W~ cross 
sections independently. 
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Rw± (stat) ± (syst) 


M < 0.5 
0.5 < \ Ve \ < 1.0 


4.3 ± 0.7 ± 0.3 
2.9 ± 0.5 ± 0.2 



TABLE VIII. (Color online) Measurements of the W cross 
section ratio, Rw, for the two e ± pseudorapidity bins. 
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FIG. 12. (Color online) W cross section ratio, Rw, for the two 
e ± pseudorapidity bins. Theory calculations at NLO from the 
fewz program using the MSTW08 and CTEQ 6.6 PDF sets 
(with 90% confidence level error eigenvector uncertainties) are 
shown for comparison. 



Our results for the measured cross section ratio are 
listed in Table |VIII| Figure [12] shows the cross section 
ratio as a function of \rj e \, where the statistical and sys- 
tematic uncertainties of the data have been added in 



quadrature. Also displayed in Fig. 12 are theoretical cal- 
culations of the cross section ratio computed with the 
fewz program at NLO. Both the MSTW08 and CTEQ 
6.6 PDF sets were used to compute the ratio; the error 
bands shown are the 90% confidence level error eigenvec- 
tor uncertainties. The predictions agree with the mea- 
sured values within the large uncertainties, which are 



dominated by the statistical precision of the W yield. 



VIII. SUMMARY 

We have presented measurements of the W + —> e + v ei 
W~ — s- er u e , and Z/-f* —> e + e~ production cross sec- 
tions in proton-proton collisions at y/s = 500 GcV by the 
STAR detector at RHIC. Theoretical predictions based 
on pQCD calculations are in good agreement with the 
measured cross sections. In addition, a first measure- 
ment of the W cross section ratio is presented. Future 
high statistics measurements of the W cross section ra- 
tio at RHIC will provide a new means of studying the 
flavor asymmetry of the antiquark sea which is comple- 
mentary to fixed-target Drell-Yan and LHC collider mea- 
surements. 
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